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The nonextreme Dp-brane solutions in type II supergravity (in the near-horizon limit) are ex- 
pected to be dual to (p + l)-dimensional noncompact supersymmetric Yang-Mills theories at finite 
temperature. We study the translationally invariant perturbations along the branes in those back- 
■ grounds and calculate quasinormal frequencies numerically. These frequencies should determine the 

f*"^ ' thermalization time scales in the dual Yang-Mills theories. 
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I As with any dualities, gauge/gravity dualities or AdS/CFT dualities (anti-deSitter/conformal field theory dualities) 
. are interesting in two respects, are interesting in two respects. On one side, supergravities give information of dual 
Q\ ' Yang-Mills theories in strong coupling regimes such as confinement. 

On another side, Yang-Mills theories give information of supergravities. Finite temperature gauge/gravity duals are 
particularly interesting in this respect. The original finite temperature gauge/gravity duality is the duality between 
Af = 4 finite temperature super- Yang-Mills theory (SYM) and type IIB string theory in the Schwarzschild-AdS5 
(SAdSs) x5 5 background Q. Thus, finite temperature gauge/gravity dualities should address long-standing puzzles 
in gravity, such as the singularity problem 0, Q (and references therein), the information paradox 0,0, and the 
Gregory-Laflamme instability 0, . 

Unfortunately, finite temperature gauge/gravity dualities have been less studied compared with the zero- 
temperature AdS/CFT. First, many evidences of finite temperature gauge/gravity dualities remain qualitative. Most 
well-known evidences are 
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I. INTRODUCTION 



1. The existence of the confinement-deconfinement transition (the Hawking-Page transition in gravity side Q) 

2. The large- N dependence of the partition function in each phase 



but quantitative understandings are still far. (The lack of supersymmetry for finite temperature is obviously the main 



obstacle.) Second, backgrounds other than SAdSs have been less studied compared with the zero-temperature cases. 1 
This paper provides one step toward these directions. The original AdS/CFT duality is motivated by the near- 
horizon limit of the extreme D3-brane. But the similar dualities are expected for the other Dp-branes (if p < 5) [Tl|. 
They are expected to dual to (p+ I)-dimensional SYM. Some qualitative features have been known for these dualities, 
but there are less quantitative studies. (Various Wilson loops have been computed, e.g., see Ref. an d references 
therein. At the zero temperature, the correlation functions are computed in Ref. |l3||.) 

Our aim is to compute the quasinormal (QN) frequencies in these backgrounds, ft is an important concept in 
black hole physics and has been widely discussed in the literature. Moreover, the QN frequencies of an AdS black 
hole have an interpretation in the dual gauge theory. Such a black hole corresponds to a thermal state in the gauge 
theory. The QN frequencies measure how perturbations of black holes decay. In the dual theory, this corresponds 
to the process where a perturbation of the thermal state decays and the system returns to the thermal equilibrium. 



'Electronic address: kmaeda@kobe-kosen.ac.jp 
^Electronic address: makoto.natsuume@kek.jp 
t Electronic address: okamura@ksc.kwansei.ac.jp 

1 For zero temperature, many backgrounds with less supersymmetry are known by adding perturbations to the Af = 4 SYM. S om e 
examples are the Klebanov-Strassler background and the Polchinski-Strassler background, which are dual to certain Af = 1 SYM 
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Thus, ON frequencies give the prediction of the thermalization time scale for the strongly-coupled gauge theory 
|14 lla lialm lialia ? It has been also argued that these modes govern the behavior of gauge theory plasmas. (See, 
e.g., Refs. |2(1 l2ll 1221 I23I l24j and references therein.) 

The plan of the present paper is as follows. First, in the next section, we review gauge/gravity dualities for Dp- 
branes. In section HTT1 we present the perturbed equations which are translationally invariant along the brane, and 
briefly review a simple numerical method 0] to obtain QN frequencies. In section IIVI we discuss the numerical 
results. We conclude in section with a summary of our results. In the Appendix, we briefly review QN modes for 
readers not having sufficient background in them. 



II. GAUGE/GRAVITY DUALITIES FOR Dp-BRANES 
A. Bulk geometry 



In this Section, we quickly review gauge/gravity dualities for Dp-branes. The relevant part of type II supergravity 
action is given by 



S = 



1 



-G 



-2<f> 



(i2 + 4(V0)' 



1 



Ft 



2(p + 2)! p+2 



The nonextreme Dp-branes are written as 

ds 2 = z-^ 2 {-hdt 2 + dx 2 p ) + zy 2 (hr x dr 2 + r 2 dn 2 _ p ), 



g 2 e~ 2 * 



(2.1) 



(2.2) 



where g s = e^°° and the harmonic functions are given by 

Mr) - < r ^-' 



r-7-P 



r p - ~ g s Nl 7 s ~P, 



(2.3) 
(2.4) 

Let us take the "decoupling" limit or the "near- horizon" limit. In order to take the limit, we assume r ~ r <C r p . 
Then, Z p = 1 + (r p /r) 7 ~P -» (r p /r) 7 "P. 
Introducing a new coordinate f, 



the metric becomes 
ds 2 



5 — p 



3 — P f ^5 — p 



-H(f)dt 2 + ( y ) dx 2 p 



~\ 2 

r 



df 2 

W) 



(2.5) 



(2.6) 



where I — r p . H and <j> are rewritten as 



H(f) 



g 2 s e- 2 * 



~\ 2 

r 



r Q \ 5 - p /r^ 
Tj (l 



(7-p)(p-3) 

5-pl\ 5 - 
2 f 



(2.7) 



It is clear that the metric is conformal to AdS p +2 x S 8 p asymptotically if p < 5. Solutions with p > 5 do not have a 
positive specific heat, so we consider p < 5. We shall use such a "AdS-frame" to calculate QN frequencies. 



2 As with all dualities, gravity and gauge theories do not have an overlapping region of validity. Our results should be regarded as the 
strong-coupling prediction of gauge theories (in the region where D-brane description is valid. See Sect. ITT^l ) 
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As is well-known, the above metric is not the SAdS solution even if p — 3. The SAC1S5 solution is given by 

ds 2 = - (ft + 1 - ^) dt 2 + + 1 - A) " + r^l (2.8) 

SAdSs has the horizon with the topology of S p , whereas the metric Ij2.fi|l has the horizon with the topology R p . We 
call such solutions "planar black holes." The p — 3 planar black hole corresponds to the large black hole limit of 
SAdSs. In order to reach the planar black hole from SAdSs, rescale the coordinates 

t — > t/a, f — ► ar, and fo — * af<j. (2-9) 

Then, the S 3 radius is proportional to a, so a 2 dil 2 ~ dx p 2 for large a. This limit, the planar SAdSs black hole, is 
invariant under the above scaling. 

We consider such planar black holes from the following reasons. First, AdS black hole solutions with S p topology 
are not known when the dilaton is nontrivial. Second, in the SYM description, SAdS corresponds to a compact SYM 
on S p , and a planar black hole corresponds to a noncompact SYM on R p . The motivation to consider a compact 
SYM in Ref. [ij is to break the scale invariance l|2.9|l . Without breaking the scale invariance, one cannot see the 
confmcmcnt/dcconfincment transition in gauge theory. Here, there is no scale invariance due to the dilaton. So, it is 
not clear whether one should consider a compact SYM. 3 



B. Validity of supergravity descriptions 

To discuss the validity of supergravity description 1|2.6I) . it is convenient to introduce SYM variables. The SYM 
coupling in terms of string variables are (See, e.g., Ref. |25j) 

g 2 YM = (2n) p - 2 g s l p -\ (2.10) 

where gyu is the (p+ l)-dimensional SYM coupling constant. The effective (dimensionless) coupling of SYM theories 
is 

(\ p— 3 
p) ■ (2- 11 ) 

The SYM perturbation theory can be trusted in the region 

9eff « 1. (2.12) 

On the other hand, one can trust supergravity solutions if both the curvature (in string metric) and the dilaton are 
small. Since 

e+~g%y )l2 lN , l 2 s R~l/g eff , (2.13) 

these conditions imply 

l«g L 2 // « A^ 4/(7 - p) . (2.14) 

Clearly, the perturbative SYM and supergravity descriptions do not overlap. For p < 3, this gives the following range 
of r (not the AdS-like coordinate f): 

(5ya/A0 1/(3 " p) ^ 4/(3 ~ p)(7 ~ p) « r/l 2 s « {g 2 YM N) 1/{3 - p) ■ (2.15) 

(For p > 3, replace the -C signs by 3> signs.) The left-hand side and the right-hand side of these inequalities come 
from the dilaton and the curvature, respectively. They have a diverging dilaton at r — and a curvature singularity 
at r = 00. One necessary condition to satisfy the above condition is N ^> 1. 



3 It is not known if there is a confinement/deconfinement transition in these theories. The dual geometry seems to suggest that there is no 
such a transition in these theories as well; namely, the specific heat is always positive for 12.61 . so there is no sign of thermal instability. 
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When the radial coordinate r is outside the region, different theories (such as M-theory) take over the type II 
descriptions. The radial coordinate has the gauge theory interpretation as the energy scale. The phase diagrams 
are discussed in Ref. For type IIA-branes, the M-brane description take over at small radius. As argued later, 
the M-brane description often reduces to a SAdS black hole, and there have been extensive works on the subject; 
this limit is relatively well-known. The other limit is the perturbative SYM description at large radius; again this 
limit is rather well-known. Therefore, we focus on the intermediate energy scale where type II supergravity is a valid 
description. 

In order to calculate QN frequencies, one places boundary conditions both at the horizon and at infinity. We 
henceforth consider the case where the horizon radius lies inside the region (|2.15|l . 4 In the large- N limit (N — » oo and 
9ym —> with a fixed large gy M N), one can enlarge this region as large as one wishes. Type II description is not a 
good description at infinity as well. One might put a boundary condition at a large radius but within the region of the 
validity. We here assume that such a boundary condition at large radius does not affect the results significantly. We 
return to the issue of the validity and discuss how different boundary conditions may affect our results in Scct. UvTI 



III. NUMERICAL APPROACHES FOR QN FREQUENCIES 
A. Basic equations 

After the conformal transformation of the metric (|2.2|l . 



ds 2 = I h —JP-\\g s e-4>)^ds 2 = Gmn dx M dx N , (:',.! ) 



the action j|2.1|) is transformed as 

s / d w x\ -a 

up to a constant, where 



e~ 



^2(p+l) -2a<t> _ 



(3.2) 



2(p + 2)! p+2 

p-3 (p-lXp-4) 5-p ^ 

a= —p> b= (7- P ) 2 ' ^ = — 9s ■ (3 - 3) 

Since we are not interested in the perturbation on S 8 ~ p sphere, we assume an ansatz for the metric: 

ds 2 = g {x p )dx»dx u + (^—^ l) dQ 2 (3.4) 



where A = — jt-^-^St^ and greek indices run from to p + 1. After the compactification, one gets a (p + 2) 
dimensional action [2|| 



S = J d p+2 x ^f~g_ e 2 ^ (R-2A + 4b (V0) 2 ) , (3.5) 



up to a constant. The p = case reduces to two-dimensional gravity coupled to a scalar, so the system locally has no 
dynamical degrees of freedom. Hereafter, we focus on 1 < p < 4. 

For simplicity, we shall only consider the perturbations which are translationally invariant along the brane. Then, 
one can set the metric as 

ds\ p+2) = g v dx"dx u = (2) g ab {x c ) dx a dx b + e 2C(xa) dx 2 , (3.6) 



4 As a matter of fact, our results are valid even for black holes with smaller horizon (for p < 3). This is because the change of a supergravity 
description (e.g., from the type IIA supergravity to the 11-dimensional supergravity) does not change the results. On the other hand, 
the horizon must be sufficiently smaller than the right-hand side of Eq. 12.151 . This is because the supergravity description must be 
valid not only at the horizon, but also in the region where QN modes decay. These issues are discussed in Sect. IWUI 
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where x a — (t,f). And the action (|3.5() becomes 



S = 



2 x ^ 



P(P~3) 2 f { 
9 



9-p 



(2 



where 



(7-p)(p-3) 



2a < 



■PC 



(3.7) 



(3.8) 



The system is two-dimensional gravity coupled to two scalars, ip\ and ip2, and there is only one dynamical degrees 
of freedom. We need to find the dynamical degrees of freedom and obtain its perturbative equation around the 
background solutions given by Eas. H2.6|) and 1)2.7(1 . Fortunately, since the background solution for ipi is constant, it 
is easy to show that its perturbation 5ipi is gauge invarint at the perturbative level [27| , and the equation of motion 
is given by 5 



= d a 



(3.9) 



where the bold face letters are background quantities. 

One can show that this equation is invariant under the scaling (|2.9|> with I fixed. By the same argument as in 
Ref. this scale invariance means that QN frequencies are proportional to fo, or the black hole temperature 

T oc fo// 2 . Hence, it is convenient to introduce dimensionless coordinates r = f o t/l 2 and t] = f/fo, and 



>2/2 



Then, Eq. (|3.9() is rewritten by 



where 



-0 2 T Z = (-dl + V ) Z = [ - Hd^Hdr,) + V ] Z , 



_ 9-p 

and the tortoise coordinate 77* is defined by 



V = vH (v+ 1 
2(7 -p) 



(1 = 



drj/H 



(3.10) 

(3.11) 

(3.12) 
(3.13) 

(3.14) 



The potential V is monotonically increasing function of 77 and V \ v= i — 0, so the potential V is positive-definite 
outside the horizon (77 > 1). 



B. A numerical method to obtain QN frequencies 

Following Horowitz and Hubney's method 0], let us calculate the QN frequencies for Dl, D2, and D4-branes. 6 
For the later convenience, we shall define ip as 

Z = e- l£jT {\ - ^-ly^/i ijj(rj) , (3.15) 

and introduce a new coordinate x = ?/~ 2 ° '> '( 5 ~p), where er = 1 for p = even and a — 2 for p = odd. When p < 5, 
the horizon and the infinity correspond to x = 1 and x — 0, respectively. Then, one gets a Fuchs-type differential 
equation: 

s(x)^(x) + Mr ^ ^(x) + V(x) = , (3.16) 

da;' i 1- 1 «i (x — iy 



5 This equation is equivalent to the massless scalar field minimally coupled to ten-dimensional Einstein metric. 

6 For p = 3, the QN freq uencies for the planar black hole correspond to the frequencies for the large black hole limit of SAdSs, which 
were calculated in Refs. Il4ll2ij . 
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x (T-p)/a_ V 2 



where 

s(x)=\x - — - j , (3.17) 

T T (7— P)/a _ 1 r -i 

t(x) = - 5-2a-p+\9 + 2cr-p-2i(5~p)Lj}x ( - 7 - p ^' T , (3.18) 

2cr x — 1 L J 

u ( x ) = - (y^J^J [ (£> + * f ) 2 x 2 < 7 - p )/ CT — (H-2*i/o>) x( 7 - p )/ <T - £ 2 a;( 5 -P)/ CT + ^ + 1)] . (3.19) 

We solve these equations by expanding ip around the horizon (x = 1). For the power expansion of ip about x — 1 
to be applicable up to the asymptotic infinity (x = 0), the radius of convergence must reach x = 0. Let us examine 
the singularity structure of Ea. H3.16|) on the complex x-plane. Equation (|3.16|l has regular singular points when x = 
and x^ 7 ~ p ^' J = 1, namely x = e 2m ' f= p (m — 1, 2, • ■ • , 7 — p) for p = even and a; = e 27Tl (m = 1, 2, • ■ • , ■^ £ ) for p = 
odd. So, the nearest singular point of x = 1 is x = and the power expansion of ip about x = 1 is applicable up to 
the asymptotic infinity (x — 0). 

The QN modes are the solutions of Ea. (|3.16|l with the following boundary conditions: 

(i) The scalar wave is purely ingoing near the horizon, Z ~ e —i<i>( T +>?»), 

(ii) The scalar wave decays at infinity, Z ~ e~ %UT rf^ l ' +1 \ (Another mode diverges.) 
The purely ingoing mode is expressed near the horizon as 

Z ~ e ~^ e -^, = e -i*T e -iu,fdr,/H _ e -^r^ _ q-iCi/q _ e -^( X _ J? -9)-«S/ 9) (3 20 ) 

which is just the same prefactor in front of ^ in Eq. I|3.15|l . So, the solution satisfying the condition (i) has the form 

00 

iP(x)=J2 a n(x-l) n , (3.21) 

n=0 

where ao is a non-zero constant. Then, the coefficients a n are obtained by the following recursion relation: 

n-l 



a,, 



n(n — 1) sq + nto 



^ [ - 1) s„_fc + fc<„_fc + u„_ fe ] afe , (3.22) 



fc=0 



where s„, i„, and u„ are the n-th order coefficients of the expansions around x = 1, e.g., s(x) = X^Li s « — l) n - 
Since the equations are linear, the coefficient ao is a free parameter (we set ao = 1). 

In order to find the QN frequencies, we need to find the solution satisfying the latter boundary condition (ii): 

00 

#)) = 5>„(-l) n = , (3.23) 

which gives a polynomial equation of uj. 

Before we solve the above equation of uj for p = 1,2, and 4 numerically, we comment on some general properties of 
QN frequencies: 

• QN frequencies are symmetrically distributed with respect to the imaginary axis of the complex d>-plane, i.e., 
there is a symmetry uj <-» — uj* . 

• The imaginary part of the QN frequency is negative so that the system in our interests is stable. 

The former statement is proved as follows: Consider any QNM, Z& with a QN frequency uj. Then, (Z&Y is also 
a solution of Ea. (|3.11|l with the frequency uj — > — uj* . Furthermore, its asymptotic form is (Z&) ~ e _l< -~" )( T +»)») 
near the horizon, and (Z&) ~ e~^~" ) r ^-(' y + 1 ) near the infinity. This means that (Z&) is also a QNM whose QN 
frequency is —uj*. 

As usual, the latter statement is proved by the "energy integral" 0, thanks to the positivity of the potential V 
outside the horizon, 77 > 1. 
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IV. DISCUSSION OF RESULTS 

To solve Eq. (|3.23() numerically, we find zeros of a partial sum i/jn = Yl n =0 a «(~l)" f° r a large N using MATHE- 
MATICA. To obtain an accurate value of QN frequencies, we need to compute on the order of TV = 100. 7 

We present the QN frequencies for p = 1, p = 2, and p = 4 cases in Table U and their distribution on the complex 
(2>-plane in Fig. ^ The numerical values of QN frequencies are normalized by the dimensionless Hawking temperature 
T = (1 1 '4ir)dH j 'd/q\ v —x — q/An. As is apparent from FigQ] the QN frequencies for each p are distributed on the 



Mode 


lor/T 


p=l 


UJl/f 


&r/T 


p = 2 

LOl/f 


Ur/T 


p = 4 


U)l/f 





7.747 




-11.158 


8.710 


-10.260 


10.488 




-5.472 


1 


13.242 




-20.594 


14.775 


-18.509 


16.232 




-8.709 


2 


18.700 




-30.023 


20.780 


-26.743 


21.805 




-11.889 


3 


24.149 




-39.450 


26.770 


-34.972 


27.319 




-15.051 



TABLE I: QN frequencies for p = 1, p = 2, and p = 4. ujr and ilii are real and imaginary parts of il>, respectively. 
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FIG. 1: QN frequencies for Dl(x), D2(o), D3(A) and D4-branes(+). Each straight line is determined by the least-square 
fitting. 

straight line. The QN frequencies are approximately given by the modes n (n = 0, 1, 2, • • • ) of the formula below: 

uj n /T = a p n + P p , (4.1) 



7 For example, for p = 2, the N = 90 results differ from the N = 100 results by about 10 5 . As a check, we have also applied the 
numerically stable continued fraction method presented by Leaver |28| and checked that the numerical values coincide with the ones 
obtained by MATHEMATICA with a good accuracy. 
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where 

{5.45 - 9.43 i for p = 1 
6.00- 8.23 i for p = 2 
6.32- 6.29 i for p = 3 ' 
5.52- 3.16 i for p = 4 

This property that the QN frequencies are approximately evenly spaced with n is numerically observed for the scalar, 
vector, and gravitational perturbations on the SAdS4 black holes [1*3] and on the SAdSs black holes [22L I23L p4| . 
The property is analytically shown for a minimally coupled massive scalar |Tflj and the vector perturbations [23| on 
the SAdSs black hole. For highly-overtone QN frequencies, one can make a comprehensive research analytically |29j|. 
However, its origin and significance are far from obvious. 



7.80- 11.17* for p = 1 

.76 - 10.27 i for p = 2 

9.88 - 8.66 i for p = 3 

10.66- 5.53 i for p = 4 



(4.2) 



A. p — 1 and p — 4 

For both p = 1 and p = 4, the dilaton gravity action l|3.5|l is simply written by 

S ~ /" d p+2 :r e 2 ^ - 2A) . (4.3) 

So, under the metric ansatz: 

4(5-p) 

one can embed the (p + 2)-dimensional action l|4.3l) into the (p + 3)-dimensional action: 

5 ~ 1 ^ +3 xy=^^ ( - 2A ) . (4.5) 

For p = 4, the embedding can be understood as the M-theory embedding (M5-brane) of the D4-brane, and this 
M5-brane reduces to the planar SAdS7 black hole (in the near- horizon limit). 

For p = 1, Dl-brane does not have a M-theory embedding because it is a type IIB object. However, various dualities 
relate the Dl-brane to the M2-brane. Under the four-dimensional pure gravity theory, the metric l|4.4|l becomes 

rf,s 2 = -(f 2 - f 3 /r)dt 2 + {r 2 - f Q 3 /r)~ l dr 2 + f 2 {dx 2 + dz 2 ), (4.6) 

where we set I = 1 for simplicity. This metric corresponds to the planar SAdS4 metric. So, one can compare our 
p = 1 results with the ones for the SAdS4 in the large black hole limit. Table [H] is the comparison with the results by 
Cardoso, et.al. [l]| corresponding to the £ = 2 gravitational perturbations (even parity) of a large SAdS4 black hole 
(1 = 1 and fo = 100). Their results are normalized with respect to the Hawking temperature T = (q/A-K) (fo/l 2 ). One 
can easily see that our results agree well with the SAdS4 results in the large black hole limit. 







Ours 


Cardoso, et.al 


Mode 


ujr/T 


til It 


ujr/T 


loi/T 





7.747 


-11.158 


7.748 


-11.157 


1 


13.242 


-20.594 


13.244 


-20.591 


2 


18.700 


-30.023 


18.703 


-30.020 


3 


24.149 


-39.450 


24.153 


-39.446 



TABLE II: QN frequencies corresponding to £ — 2 gravitational perturbations (even parity) of a large Schwarzschild-AdS black 
hole with fo = 100 and 1 = 1. 



On the other hand, the p — 2 case is not related to a higher-dimensional SAdS, contrary to the p — 1 and p = 4 
cases. Obviously, one can always embed type IIA objects into M-theory. However, the resulting geometries are not 
SAdS black holes, and QN frequencies for such geometries are unknown. 

The D2-brane is embedded as a M2-brane, but the embedding of the geometry (I2.6|) is not the SAdS4. This is 
because the embedding corresponds not to the standard M2-brane, but rather corresponds to the so-called "smeared 
M2-brane." In a sense, we have obtained QN frequencies of the "smeared M2-brane" via the D2-brane. 

Even the smeared M2-brane description is not valid at lower energy. The smeared M2-brane becomes unstable at 
lower energy due to the Gregory-Lafiamme instability and decays into the M2-brane on a circle. Then, the SAdS4 
calculation suffices for such a small black hole flfJI . 



C. Sensibility on the boundary condition 

We placed the Dirichlet condition at infinity to compute QN frequencies. However, supergravity description often 
breaks down at infinity. So, strictly speaking, one must place an ultraviolet cutoff and put a boundary condition at 
the large finite radius. Here, we discuss how such a boundary condition may change our results. 

First of all, it is not clear what boundary condition one must impose. The out-going wave certainly leaks out to the 
asymptotic infinity. So, a simple Dirichlet condition at the cutoff does not suffice. However, appropriate boundary 
condition is not clear, so here we use a Dirichlet boundary condition for illustration to see if our results are sensitive 
to the boundary condition. 

FigureElis the result of QN frequencies (for p — I) by imposing the Dirichlet condition at various radii Xb = (?"u /rb) 2 . 
(Here, r& is the location of the boundary condition in the original D-brane coordinate r.) As the figure shows, the 
result begins to converge for a sufficiently large boundary radius (say, Xb < 0.02 or rb > ro/\/0.02). For the boundary 
condition placed in the plateau region, the result is effectively the same as the one with the boundary condition at 
infinity. In other words, the result is insensitive to the boundary condition. 




0.02 0.04 0.06 0.C 



0.1 0.12 



0.1 0.12 



FIG. 2: The dependence of the lowest QN frequency on the location of the boundary condition Xb- The left figure (right figure) 
is the real (imaginary) part of the QN frequency. 

Let us reinterpret the result in terms of SYM variables. The Dl description is valid for gynN 1 ^ 6 <C r/l 2 <C 
((/y M A^) 1 / 2 in SYM variables. Actually, gravity description is valid even inside the infrared cutoff. The type IIB 
fundamental string description takes over (via type IIB S-duality), but the gravity computation is the same as the 
DI case. Supergravity descriptions are valid as long as gym <C r/l 2 <C (gy^jN) 1 ^ 2 . 

The theory has two independent parameters gyu and N. The infrared cutoff can be controlled by gym whereas 
the ultraviolet cutoff can be controlled by N. Thus, the condition Xb < 0.02 must be satisfied for a large enough N. 
This is indeed true. 

In order for the QN frequency to be insensitive to the boundary condition, must lie within this region. This 
implies g YM < r /l 2 < (0.02g YM N) 1/2 . This condition has a solution if N > 50. 8 



If one chooses a boundary condition different from the Dirichlet boundary condition, the conclusion could be totally changed (even 
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V. SUMMARY 

We have computed the scalar QN frequencies for asymptotically AdS black holes (in appropriate frame), which 
correspond to the decoupling limit of nonextreme Dp-branes (for p < 5). We consider the translationally invariant 
perturbations along the branes. The dual gauge theory is described by (p+ l)-dimensional super- Yang-Mills theories 
at finite temperature. 

As discussed in Sect. I1VI one can embed the system into higher dimensions for p = 1 and p = 4. Then, the bulk 
geometries become just SAdS black holes. Thus, these cases can be reproduced by QN frequency calculations of 
standard SAdS black holes. (However, they correspond to different region of the validity, so one must be careful to 
its interpretation. For example, the p — 4 case reduces to the M5-brane, but this must be interpreted as the recovery 
of conformality at high energy.) This gives a nice check of our approach since our perturbation is more involved 
compared with the standard cases. (Only the minimally-coupled test scalar field is often considered to calculate QN 
frequencies.) We explicitly checked that our results coincide with the SAdS results. 

On the other hand, p — 2 case is not obtained from higher dimensional SAdS black hole. Thus, the embedding of 
the p = 2 case does not help to simplify the calculation. 9 The fact that QN frequencies are evenly spaced had been 
observed for pure gravity black holes. Our p = 2 result may indicate that this is true even for some dilatonic AdS 
black holes. 

It is difficult to calculate the thermalization time scale in the gauge theory. But Ref. pointed out that the time 
scale is likely to be independent of the 't Hooft coupling. A free field theory never thermalizes and at weak coupling 
the time scale should be very long. So, clearly the time scale depends on the coupling at weak coupling. However, in 
the strong coupling, the relaxation time scale is expected to be the order of the thermal wavelength. 

In the dual supergravity description, this means that the QN frequencies depend only on the Hawking temperature. 
This expectation is indeed true for SAdS black holes [14(, and we found that it is also true for Dp-branes. (Table [I] 
shows that QN frequencies are linear in temperature.) 

Finite temperature gauge/gravity dualities deserve further study. To qualitatively check the duality, one wishes 
to check the time scale in the dual theory. Currently, it is difficult to calculate the time scale in gauge theories. 
Such a calculation has not been carried out even in the SAdS cases. But gauge theory understanding is essential to 
solve the long-standing puzzles such as the singularity problem, the information paradox, and the Gregory-Laflamme 
instability. 

Note added: After this work was completed, we were informed about work by Iizuka, Kabat, Lifschytz, and Lowe 
|30|. where they also computed QN frequencies for nonextreme Dp-branes. The differences from our paper are as 
follows: (1) We computed not only the lowest modes, but also higher modes as well. As a result, we are able to 
see how eigenvalues are distributed in the complex plane. (2) The perturbation considered in the paper eventually 
becomes the same as ours. However, they regarded the perturbation simply as a minimally-coupled test scalar field, 
whereas we derived the equation from a combination of the gravitational and dilaton perturbations, which is part of 
type II spectrum in those backgrounds. (3) We also consider the validity of supergravity descriptions. 

We would like to thank D. Kabat for communicating their results. 
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APPENDIX: REVIEW OF QUASINORMAL MODES 

We briefly review quasinormal modes (for a comprehensive review, see |3ljh 



the existence of the plateau region). So, one should not take the conclusion seriously. Our point here is that there exists at least one 
boundary condition, where the Dirichlet condition at infinity gives a sufficient accuracy. 

This fact also applies to the p = case. The DO-brane is a M-theory Kaluza-Klein state, so this corresponds to the near-horizon limit 
of a pp-wave. 
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Let us consider the initial-value problem of a linear wave equation for a scalar field $ 

(-n + V(x) )$(x) =0 . (A.l) 
As is well-known, we can formally solve the initial-value problem using the retarded Green function G R as [3^ | 

$(x) = / dE(«) [ (n^G^z)) $(*) (n"fl£*(*)) ] , (A.2) 



where dS(z) is the infinitesimal surface element of the initial surface So and is future-directed unit normal vector 
to S - 

The retarded Green function G R (x ; z) in [d + l)-dimensional spacetime is the solution of the inhomogeneous wave 
equation 



( - D x + V(x) ) G R (x ; z) = -5^ d+l \x - z)/^gjx) , (A.3) 

satisfying the causal condition, G R (x ; z) = for x J + {z). 

For static spacetimes, ds 2 = N 2 (x) ( — dt 2 + 7^ {x)dx i dx^ ) , so the Green function is time-translationally invariant, 
Gr(x ; z) —: G R (t, x ; r, z) — G R (t — t, x ; 0, z), and it is convenient to use the "frequency-domain Green function" 
gn defined by 



/>oo 

g R (x;z\s) = [N(x) N{z)] (d - l)/2 / dt e~ st G R (t, x ; 0, z) , 

Jo- 



(A.4) 



which is simply the Laplace transform of Gr (with the weight N^ 1 ^ 2 inserted for convenience). If there exists such 
a Laplace transform g R , an abscissa of convergence a exists and g R is well-defined for 5R(s) > cr. Then, g R satisfies 
the equation, 



[s 2 -A x + -- - ] g R (x;z\ S ) =-&<■*> (x - z) / yffix) , (A.5) 

where A is the Laplacian with respect to the metric 7y and the dots denote the terms by the "effective potential" . 
After solving Ea. ljA.5p with a suitable boundary condition, G R is obtained by 

/><x+i 00 1 

G R (t,x;0,z) = [N(x) N(z)f- d)/2 / ^ e st g R (x;z\s) . (A.6) 



One can define g R for 3?(s) < cr by the analytical continuation from SR(s) > cr. Provided that |s e s * <7#| — — for 
5ff(s) < a, the contribution to G R in Ea. (|A.6|l comes from the singularities of g R in 3?(s) < cr. The pole singularities 
of <?,r, {s„}„=o,i,2,-, contribute to G R as 

G R (t,x;r^) = [iV( a; )iV(z)]^/ 2 ]T e s " Res[ g H (x ; z\s) ; s„ ] + • • • , (A.7) 

n 

and the corresponding frequency for each pole u> n = is n (~s{u) n ) < cr) is called a quasinormal frequency. So, a 
quasinormal frequency is given by a pole of the retarded Green function in the frequency domain, g R (x ;z\s). 

In order to obtain g R , one must pay attention to the boundary conditions for Ea.l |A.5| l. We consider two cases 
separately. 

(i) For asymptotically flat black hole spacetimes, there are two asymptotic regions, the near-horizon region and 
spatial infinity. Equation l|A.5fl often has simple forms in these regions: 



(s 2 -dl)g R {x;z\s)~0, (A.8) 

where we consider a massless scalar field for simplicity. The tortoise coordinate of a;, r*, is defined such that 
the spatial infinity (the horizon) corresponds to — * 00 (—00). There are two independent solutions e ±sr " in 
each asymptotic region. 

Recalling that g R for 3?(s) < cr is given by the analytic continuation from 9?(s) > cr, we first consider the 
boundary conditions of g R for 5R(s) > a. Since diverging solutions are unphysical, g R for 3?(s) > a should 
behave as 

g R (x;z\s) — > e sr * (e~ sr *) for r* -> -00 (00), (A. 9) 
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so that 

G R {t,x;0,z)~e st ~g R {x;z\s)^e s{t+r * ] (e^*^^) for r, -oo (oo) . (A.10) 

This means that the appropriate boundary conditions for the retarded Green function G R in this case are a 
purely ingoing wave near the horizon and an outgoing wave at the infinity. Since g R for 3?(s) < o is defined by 
the analytical continuation from 5ft(s) > a, g R for 5ft(s) < a also satisfies the same boundary conditions (|A.9(I . 
in spite of the diverging behavior. 

(ii) For asymptotically AdS black hole spacetimes, the tortoise coordinate 77* is defined such that the null infinity 
(the horizon) corresponds to 77* — > (—00). Although Ea. lA.5|) has the same form as Ea. l|A.8|l near the horizon, 
Eq. (|A.5|) has a different form near the null infinity, 

(-81 + C/r) 2 )g R (x;z\s) ~0 for 77* -» , (A.ll) 



1/2 I a/ ^*~|~ 1/4 

where C is a positive constant. Thus, for 77* ~ 0, we obtain y^(a;;z|s) ~ 77* or G R (t, x ; 0, 2) ~ 

^rf/2± v / c+i/4 _ : ^^±^ wnere we use N^^ d ^ 2 (x) oc 77» d 1 '^ 2 for 77* ~ 0. Since we have A_ < < A+ in many 
cases, we henceforth assume A_ < < A + . 

Again, we first consider g R for 3?(s) > a. Since diverging solutions are physically unacceptable, g R for 3?(s) > a 
should behave as 

{& s ^* for tj^ — ^ — 00 

i/2+VcTT7i . ( A - 12 ) 
77* tor 77* — > 

so that 

st a + f V * n ■ (A.13) 
e 77* for 77* — > 

We should set the same boundary conditions HA.12|) to g R for 3?(s) < cr because 5^ for 3?(s) < cr is defined by 
the analytical continuation from 5R(s) > a. 

To summarize, quasinormal frequencies are obtained by finding the poles of g R satisfying Eq. I|A.5|) with the boundary 
conditions. However, it is easy to show that, for a pole of g R , s n , there exists an eigenmode 4> n (x) of the homogeneous 
equation of Eq. (|A.5|I 

(s 2 - A + ---) 4> n (x) = , (A.14) 

which satisfies the boundary conditions (|A.9|I for asymptotically flat black hole spacetimes 31], and l|A.12() for asymp- 
totically AdS black hole spacetimes. Thus, in practice, one can also find the quasinormal frequencies by solving this 
eigenvalue problem. In the main text, we solve this eigenvalue problem (|A.14|) with the boundary conditions l|A.12|) . 
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